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Abstract 



We derive a symmetry property for the Fourier-transform of the colhsion- 
C/^ ■ less sound modes of Bose condensates in anisotropic traps connected with a 

■ somewhat hidden conservation law. We discuss its possible observation by 

I dispersive light scattering. 
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O I low temperatures numerous aspects of these systems have been investigated both from the 

^ ' experimental and the theoretical side. Adding a time dependent perturbation to the trap 



Since the achievement of Bose-Einstein condensates in trapped alkali gases ||T]-pi at ultra- 



potential the lowest frequencies of the density oscillation spectrum were measured by observ- 
ing condensate shape oscillations either by time of flight measurements or by dispersive 
light scattering 0. In this paper we want to investigate these mode functions in detail and 
I show how they reflect a special symmetry of the system. We disuss how this symmetry 
property could be observed in phase-contrast imaging measurements and comment on 
O ■ the feasibility of its observation in inelastic light sacttering. 

We restrict ourselves to the density fluctuations of energies much smaller than the chem- 
^ ! ical potential {hu <C /x) of a trapped Bose condensate at temperature T = 0. The equations 
^ I governing these collisionless sound modes in a weakly interacting system in mean field ap- 
proximation can be written in the form of hydrodynamic equations P,[TO|. For isotropic P|, 



axially symmetric |T^ and even fully anisotropic harmonic traps |T^ the solutions of these 



equations can be classified in terms of three quantum numbers. We consider these mode 



functions in Fourier space where we scale all momenta by ki = ki \Jmuj'^/2^ for i = x,y,z 
with chemical potential fi and trap frequency Ui in direction i. 

Before going into any technical details we state the central result of this paper and discuss 
how it may be tested experimentally: We show in this paper that the Fourier transform 0(k) 
of the mode functions for traps of arbitrary anisotropy factorizes into a radial part depending 
only on the modulus of the scaled wave vector 



k=\k\ = ^2^/m(fc2/^2 + ^2/^2 + ^2/^2) (0.1) 

times an angular part g{6k, (pk) depending only on the two angles 6^, (pk of spherical coor- 
dinates in scaled Fourier space. Furthermore the radial part does not depend on the trap 
anisotropy anymore and can be determined explicitly. We simply have 
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0(k) = ^ X (0.2) 

where ji, denotes the spherical Bessel-function, the positive integer degree z/ depending on 
the three quantum numbers of the mode functions, which will be specified below. This 
factorization is remarkable since neither the original problem nor the problem with rescaled 
momenta are spherically symmetric. In fact it follows from a somewhat hidden conservation 
law of an operator with spherical symmetry after rescaling as we shall show below. 

Before we derive (p.2|) we wish to illustrate this property by discussing how it could be 
observed in an experiment. One possibility to observe density fluctuations is by the direct 
observation of a condensate by phase-contrast imaging A beam of light off resonace 

with any atomic transition accumulates a phase shift on its way through the condensate, 
which is proportional to the density of atoms integrated along its optical path. In phase- 
contrast imaging the direct beam of light is shifted by 7r/2 by passing it through a A/4-plate 
in the focal plane. In the interference pattern between the direct beam and the deflected 
light on a screen phase differences are thereby turned into amplitude differences. Since 
typical angles of deflection in the experiment of ^ are only of the order of several mrad 
we can approximate the paths of deflected light also as straight lines parallel to the optical 
axis. The observed amplitude differences are therefore proportional to the projection of 
the three dimensional number density of the condensate on a two dimensional screen. By 
this technique it is possible to take time resolved sequences of pictures of an oscillating 
condensate 0, after its excitation via a modulation of the trap with a certain frequency. 

By Fourier transformimg these sequences of pictures in time with the known mode- 
frequency one gets the projection of the density fluctuation mode 0(r) onto the screen, i.e. 
the mode function integrated along the direction of the optical axis. The Fourier transform 
0(k) of the mode function for any wave vector k perpendicular to the optical axis can 
then be obtained by taking the spatial Fourier transformation of this picture with the wave 
vector k in the image-plane. To test relation p.2|) the Fourier transform must be taken for 



a range of fc-values with fixed direction of k, and thus also of the scaled k, calculating the 
modulus k from ( p.lj ). The chemical potential neccessary for the application of ( |0.1| ) must 
be determined from the number of atoms in the trap. For fixed direction of k g{6k, (pk) only 
occurs as a constant factor, which of course depends on the direction chosen and which may 
vanish in some directions. 

The scaling introduced for the momenta corresponds to scaling all coordinates with the 
Thomas- Fermi radius in this direction according to = fj ^2ix/mujf for i = x,y, z. By 
this spatial rescaling an anisotropic harmonic trap is transformed to an isotropic one and 
pictures of elliptical condensates are transformed to spherical ones. Therefore, instead of 
calculating the spatial Fourier transform of pictures of mode functions in scaled momenta 
one can also first unfold these pictures to spherical symmetry and then apply normal Fourier 
transformation. In an experiment the optical axis probably would be oriented along one axis 
of the trap, e.g. the y-axis. The simplest case then certainly is to take the orientation of 



k either along the z-axis, with k = kzy2fi/mujl, or along the x-axis with k = kxy2fi/mu'^. 

For an arbitrary, but fixed orientation of k we have to keep the ratio k^/k^ constant with 
the modulus given by (|0.1|). In summary, once having taken a sequence of phase-contrast 



images for a particular mode, relation (|0?^ ) could be tested easily by a simple analysis. 
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Now let us derive the property ( p.2|) . For a Bose gas below the Bose-Einstein transition 
the usual separation of the particle field operator into a classical part 0o(r), the condensate 
wave function, and a residual operator 0(r), describing excitations out of the condensate, 
yields for the number-density operator 

p(r) = |0o(r)P + (0o(r)0(r) + 0o(r)0^(r)) + 0^(r)0(r) . 

Thermal averaging of the last term gives the condensate depletion at T = plus the density 
of thermally excited atoms. At low temperatures for a weakly interacting system the con- 
densate density |0o(r)P contains most of the particles, and the second term dominates by 
far the last term. Restricting ourselves to the low-lying density fiuctuations with huj <^ /i 
P], one single mode of these dominant density fiuctuations with mode frequency uj takes the 



form |0|,|TT 



/.*(r)0(r) 




Here g is the atom-atom interaction constant Anh'^a/m, with the s-wave scattering length 
a, and 0(r) a normalized density fiuctuation mode. These density fiuctuations and the 
corresponding velocity fiuctuations can be described by hydrodynamic equations of the Euler 
type. From this description one can derive a wave equation PJTO[| 



-c^2 0(r) = V-^|0o(r)pV0(r) (0.3) 



m 



which both the normalized density modes and the velocity fiuctuation modes have to fulfill 



Features of Bose condensates are most readily observable in the limit of large condensates. 
In this limit the Thomas-Fermi approximation for the condensate wave function 5f|(/)o(r)p = 
{jJi — U {r))Q{fi — f/(r)) is justified [Q, which restricts solutions of ( |(J.3| ) to the interior of the 
Thomas- Fermi surface fi = U{r). For condensates in isotropic harmonic traps the spectrum 
and the eigenfunctions of (|0.3|) were calculated by Stringari 0. They depend on a radial 
and an angular quantum number and /, the spectrum being degenerate with respect to 
the azimuthal quantum number m/. 

The atom-traps used in the different realisations of Bose condensates in alkali atoms 
can be decribed in very good approximation by an anisotropic harmonic potential U{r) = 
miuj'l.x^+ujyy'^+ujlz^) /2. Here we do not restrict ourselves to any axial or spherical symmetry. 
Due to the anisotropy of the potential new features show up: the azimuthal quantum number 
nil is no longer a good quantum number, at energies of the order of /x the system ceases to 
be strictly integrable, and a large chaotic component can be seen in the classical phase space 
T3| . However in the limit of low energies {huj <^ /i) the systems reduces to an integrable 



one 



16 1 and an additional constant of motion can be identified [14 



^ = -E:!^(7^)' + (r-V)^ + 3r.V. (0.4) 



muf dr. 

We shall now show how ( p.2|) follows from this conserved quantity. 
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Let us introduce a scaled position-vector r with dimensionless components fj all in the 
range [0, 1] by the scaling = fj \j2iji/'mujf for i = x,y, z . Any harmonic trap potential 
now becomes U{t) = /xf^ and the Thomas-Fermi surface is mapped to the unit sphere. The 
wave equation (p.3|) takes the form 



(f) 



Now let us change to scaled Fourier space. Since terms from partial integration don't occur 
due to the vanishing of the condensate wave function at the condensate boundary we simply 
can replace id/dkj and d/dfj — > ikj. The modulus of the scaled momentum k is given 
by ( p.l|) . After simple rearrangements we get from 



cu^ 0(k) 



E[k'fc.^(l + V,^)+..f^.i] 0(k). 



Next we introduce spherical coordinates in scaled Fourier space in the usual way as kx = 
k sin 6k cos (pk ,ky = k sin 6k sin 0^ and kz = k cos 6k . The explicit form of the wave equation 



then becomes 

,2 
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uj'^ (f){k, 6k, (pk) = I (^(t^x. sin^ 6k cos^ 0^ + sin^ 6k sin^ 0^ + ^1 cos^ 0^) 



-|- sin 26k i^x cos 0^ -|- a;^ sin (pk — lo ' 
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sm-i(()k{ujl - ujfj-—j (f){k,6k,(pk) 



We see that the radial part, depending on k only, can be separated from this equation 
which justifies the factorized form of the separation ansatz (p.2| ). Our point now is that the 
separated equation for the radial part simply is the eigenvalue equation for the operator B 
in the scaled Fourier representation. To see this let us transform B to the k-representation. 
Starting in coordiante space again, scaling the coordinates in ( p.4|) first, and then changing 
to Fourier representation we get 



B 



[f • V)2 + 3f • V 



dk dk 



(0.5) 



which is indeed precisely the operator occuring in the radial part of the wave equation. 
Remarkably the eigenvalue-equation of B does not depend on the frequencies Ui anymore, 
its spectrum and its solutions are the same for arbitrary anisotropies. 

The eigenfunctions of B in the k-representation can be found in terms of spherical Bessel- 
f unctions ji, 



B 



k 



3u{k) 
k 



The spectrum of B was previously derived in |T^, however, without stating the eigenfunc- 
tions. In the notation of [|14| for axially symmetric traps one has v — 1 = n + irti , where mi 
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denotes the azimuthal quantum number. Previous solutions of (|0.3| ) started with a polyno- 
mial ansatz in Cartesian coordinates [|,15,1^ for 0(r) . We note that z/ — 1 is simply the 
highest degree of such a polynomial solution. For the general anisotropic case the depen- 



dence of V on the three quantum numbers is given in |T^. In terms of the quantum numbers 
of the isotropic problem one has v — \ = 2nr + 1 . For this isotropic case ( p.2| ) was already 
derived by Fetter [|l^ by direct Fourier transformation of the spatial mode functions. Since 
the conserved quantity B does not depend on the trap frequencies u!x,^^y,^z, subspaces 
of the Hilbert space for different values of the quantum number u remain orthogonal for 
all anisotropies. So the same dependence on the modulus of the momentum k as for the 
isotropic case holds for all anisotropies. 

Let us remark now on the range of validity of our discussion which is limited by our 
use of the approximate low-frequency form ( p.3| ) of the wave equation and the Thomas- 
Fermi approximation. For the description by the wave equation (p.3|) we have to require 
hu fi. A second condition arises from the finite width of the boundary layer /j = 
(2jj/mujfY^'^{huJi/4:jj,Y^^ for the directions i = x,y,z, limiting the momenta k for which the 
Thomas- Fermi approximation can still be used: There is a maximum momentum ku = 27r/Zj, 
or in scaled variables ku = 27r(4/i/^j)^/^, for which ( p.2|) is valid. For the cigar shaped 
condensate of this condition is in radial direction {kir ~ 140) more limiting than in 
axial direction [kia ~ 770). Both limiting values place no restrictions on the experimental 
observation for the low lying levels, however for higher lying levels they may be relevant. Let 
us consider the case of axial symmetry somewhat more generally: For a level with quantum 



numbers (n,j,mi) in the notation of relation (|0.2| ) can be observed best in the range 
of the first maximum /cmax of the Bessel function jn+m+i{k), scaling as fcmax ~ (n + mi) for 
large n,mi. If we have n ^ mi, then uo n and both limiting conditions on the energy 
^njm ^ 1^ and on the momentum fcmax ^ h roughly coincide. However if n -C m;, we have 

1/2 

uj m/ and the condition kmax ^ h is the more limiting one. 



As an alternative experiment to determine the collective excitations and to observe ( p.2|) 
one might also think of inelastic scattering of light [|10],O. Unfortunately, this turns out 



to be hardly feasible, as we now show. Light scattered from an atomic gas off resonance 
with any atomic transition couples to the number density of atoms. The cross section of 
elastic scattering of light is proportional to the square of the spatial Fourier-transform of the 
equilibrium density |0o(r)p. For an inhomogenous system elastic scattering occurs in a finite 
angle, fixed by the condensate size. Density fluctuations can be seen as inelastic scattering 
of light. The cross section of light shifted in frequency by u via scattering from the density 
fluctuation is again proportional to the Fourier-transform of this mode function |jTO|. One 
might hope to use this to test relation ( |0.2| ). For scattering transfering the wave vector k 
and at T = we get as a ratio of elastic to inelastic cross section {N fi/hu) x |po(k)/0(k)p, 
which unfortunately turns out as particularly unfavorable both in the limit hu <^ fi and in 
the limit of large particle numbers. 

The larger the transferred momentum k, the better inelastic scattering can be distin- 
guished from the background of elastic scattering. The maximum momentum, for which 
( p.2|) is valid, is the momentum ki determined by the boundary width. For present experi- 
ments 1^,0 this is by a factor 2 — 5 smaller than the largest momentum transferable by light 
scattering, occuring in the case of back-scattering. Therefore the inelastic cross section from 
( p. 21 ) describes scattering for not too large angles of deflection only. For the maximum mo- 
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mentum permissible by our approximations the elastic cross section is of the order of l//cf, 
which follows from the (on this scale) discontinuous first derivative of |0oP at the boundary. 
For small quantum numbers the inelastic cross section from (p.2|) scales as 1/fcf. So col- 
lecting the prefactors we get as a minimum overall ratio from elastic to inelastic scattering 
(/i/^a;)(^j/4/i)^/^ for inelastic scattering from level uj in direction i. For the experiments 
of this ratio exceeds 10^ for all directions for the low lying levels. This unfavorable ratio 
probably precludes the possibility to test relation ( p.2| ) by inelastic light scattering. 

In conclusion we have shown that the scaled Fourier-transform of the hydrodynamic 
mode functions factorizes in a radial part and an angular part, even for harmonic traps 
of arbitrary anisotropy. The radial part is fixed by only one quantum number and has the 
simple explicit form ( p.2| ) . This result was derived from a previously found conservation law. 
The eigenvalues and eigenfunctions of the conserved operator B determine completely the 
radial mode function in scaled Fourier space . We have discussed the range of validity of this 
property and have suggested how it could be observed in experiments using phase-contrast 
imaging. 

M.F. thanks H. Schomerus for a valuable discussion. This work has been supported by 
the Deutsche Forschungsgemeinschaft through the Sonderforschungsbereich 237 "Unordnung 
und groBe Fluktuationen" . 
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